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Abstract 

> . 

CN . The latest g^ — 2 measurement by Brookhaven confirms the earlier measurement 

^— l ' 

■ with twice the precision. However, interpretation of the result requires specific 

assumptions regarding the errors in the hadronic light by light (LbL) correction 
and in the hadronic vacuum polarization correction. Under the assumption that 
the analysis on LbL correction of Knecht and Nyffeler and the revised analysis of 
£pH, Hayakawa and Kinoshita are valid the new BNL result implies a deviation between 

experiment and the standard model of 1.6a — 2.6a depending on the estimate of 
i-C the hadronic vacuum polarization correction. We revisit the — 2 constraint for 

mSUGRA and its implications for the direct detection of sparticles at colliders and 
for the search for super symmetric dark matter in view of the new evaluation. 



1 Introduction 

The BNL — 2 Collaboration has announced]]]] an improved result for a M = (g M — 2)/2 
with twice the precision of their 2001 analysis ||. Here we interpret these results in 
the context of supersymmetry. It has been recognized for some time that the muon 
anomalous moment can act as an important probe of physics beyond the standard model 
especially of supersymmetry. This is so because the new physics contribution to the 
leptonic anomalous magnetic moment a\ scales as mf/A 2 and is proportional to the square 
of the lepton mass. Thus provides a more sensitive probe of new physics than a e even 
though a e is more sensitively determined while the a T determination is far less sensitive 
to be competitive. Because of the above any improvement in the determination of a M has 
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important implications for new physics beyond the standard model. Thus the improved 
measurement of BNL experiment in 2001 1| which gave a e * p = 11659203(15) x 10~ 10 
resulted in a large amount of theoretical activity to explore the implications of the results. 
At the time of the BNL2001 result the standard model prediction consisting of the qed, 
electro-weak and hadronic corrections was estimated to be a^ M = 11659159.7(6.7) x 
10~ 10 [[j|] which gave a e * p — a^ M = 43(16) x 10~ 10 implying a 2.6a difference between 
experiment and theory. The above difference was based on a light by light (LbL) hadronic 
correction of 0, |5| a^ ad (LbL) = —8.5(2.5) x 10~ 10 which was later found to be in error|6|, 
Thus the analysis of Knecht et.al.g gives § a h ^ ad (LbL) = 8.3(1.2) x lO" 10 while 
the revised analysis of Hayakawa and Kinoshita gives0 a^ ad (LbL) = 8.9(1.5) x 10 -10 . A 
partial analysis of a^ ad (LbL) has also been given in Refs.[[5], [|. Thus the partial analysis 
of Ref.g finds a h ^ d (LbL) = 8.3(3.2) x 10~ 10 and the work of Ref.@ which computed the 
pion pole part finds a^ ad (LbL : ir°pole) = 5.6 x 10~ 10 . Again the sign of these corrections 
agree with the sign of the reevaluation of this quantity in Refs.|| [?[]. Using the average 
of the first two|6|, |7| which are the more complete calculations and using the hadronic 
vacuum polarization correction of Ref.|]T0| one finds a^ M = 11659176.8(6.7) x 10~ 10 and 
the BNL2001 result gives a e * p — a^ M = 26(16) x 10~ 10 . This difference corresponds only 
to a 1.6cr deviation between theory and experiment. 

Before proceeding further we wish to discuss a bit further the issue of errors in the 
hadronic corrections specifically the LbL correction and the a 2 vacuum polarization cor- 
rection. Regarding the LbL hadronic correction, in addition to the analyses mentioned 
above, i.e., Refs.|| 0, ||, [| there is also the analysis of Ref.|llJ] based on chiral perturba- 
tion theory. This analysis finds a h ^ d (LbL) = (1.3^+3.1(7) x 10~ 10 where C is a correction 
arising from the sub leading contributions. These contributions are either of 0(a 3 p 2 /A 2 ), 
where p is a mass of order m M or m n and A is a hadronic scale ~ 1 GeV, which are not 
enhanced by a factor of Nc (the number of quark colors) or of 0{Nco?p 2 /A 2 ) but are 
not enhanced by large logarithms. The coefficients C in this analysis is an unknown low 
energy constant (LEC) which contains the contributions of nonperturbative physics at 
short distance. The authors of Ref.]TT| view the parameter C as basically unconstrained 
except through the measurement of the anomalous moment itself. In our analysis we 
assume the validity of the analyses of Refs. |6], . Aside from the LbL contribution, the 
other source of error in the hadronic correction is the contribution from the a 2 hadronic 
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vacuum polarization correction ||T0|, [T2| , [TJ^ |IJ]] and this error is an important component 
in extracting the deviation between experiment and the standard model. 
The new BNL result gives the world average on a M as 0. 



a ex P = (11659203)(8) x 10 



-10 



(1) 



With the LbL correction of Refs. || [7| and assuming the leading hadronic correction to lie 
in the range 692(6) x KT 10 to 702(8) x 10~ 10 
lies in the range a™ = 11659177(7) x 10~ 10 to 
to 



10, 12, 13 1 the standard model prediction 
SM = 11659186(8) x 10~ 10 which leads 



n ex P — n SM 



- (26)(10) x 10~ 10 
to (17)(11) x 10- 10 



(2) 



which corresponds to a difference between experiment and theory of about 1.6cr to 2.6cr. 



2 Interpreting the New BNL Result for Supersym- 
metry 

Next we discuss the implications of these results for supersymmetry. The supersymmetric 
electroweak contribution to a M arises from the chargino and neutralino exchange correc- 
tions. For the CP conserving case the chargino (x ) an d the sneutrino (z>) exchange 



contributions are typically the larger contributions and here one finds fll5 
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where Al(Ar) are the left (right) chiral amplitudes and F 2 are form factors as defined 
in Ref . |TJ| . The accuracy of the result of Ref . [T5| was further tested by taking the super- 
symmetric limit of the results of Ref. [jHJ in Ref. [Tj| where in addition the effects of CP vio- 
lation on — 2 were also investigated. The amplitude of Eq. (3) is dominated by the chiral 
interference term proportional to AlAr. Since Al ~ 1/ cos/? (tan/5 =< H 2 > / < Hi > 
where < H 2 > gives mass to the up quark and < Hi > gives mass to the down quark 
and the lepton) one finds that a^ USY ~ tan/3 for large tan/3[T7l ITS]. Further, one finds 
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that Al depends on the sign of /xm2, where m 2 is 577(2) gaugino mass and [i is the higgs 
mixing parameter both taken at the electroweak scale, and thus the sign of a^ USY is de- 
termined by the sign of The above leads one to the result that over a large part of 



the parameter space unless tan/? ~ 1 one finds [If, IS 



a susY > q^ 2/i > q. 



a SUSY < 0? m 2/ u < 



(4) 



where we have used the sign convention on fx of Ref . [|1| . These results have been confirmed 
numerically in a wide class of models |Tj], [Tj| 20]. The implications of the BNL result 
of 2001 has been analyzed extensively in the literature pl|, p2; , |2^ 
framework of mSUGRA |26|[ as well as in a variety of other models. 



2~&i within the 



We analyze now the implications of the new result within mSUGRA[26j. We note in 
passing that for a class of models with large extra dimensions the correction to g-2 from 
Kaluza-Klein modes is rather small [T27I] once the constraints arising from Kaluza-Klein 



corrections to Gp are taken account of|p7|. Thus the extra dimensions do not interfere 
in the extraction of implications of the new result for supersymmetry. In our analysis we 
will use the 2.6a deviation in Eq.(2) as the default value for the deviation but we will take 
a 2a error corridor in the analysis. This error corridor includes the case corresponding to 
the smaller value of the deviation in Eq.(2), i.e., l.Qa and taking a la error corridor in 
the analysis. Our assumption thus gives us the constraint 
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< 46 x 10 
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Attributing the entire difference a e * p — to supersymmetry and using Eqs.(2) and (4) we 
find that the sign of /i is positive confirming a similar result arrived at in analyses based on 

3|. The sign of /i is of great importance for dark matter |23|, B3j] 
iSL PI. 



the BNL2001 datafP 
and for b 



t unification H 



As is well known b — > s + 7 imposes an important 



constraint on the parameter space of supersymmetric models|2^, p0 |. The constraint 
is very stringent for negative /1 eliminating most of the parameter space of models. For 
positive n the constraint is not very strong for small tan (3 but becomes a strong constraint 
for large values of tan f3. The imposition of this constraint involves the standard model 
value and the experimental value both of which have significant errors. The standard 
model prediction of this decay including the next to the leading order correction has been 
given by several authors |31|. Further, there are several recent experimental determinations 
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of this decay p2| . In our analysis we take a 2a error corridor between experiment and the 
prediction of the standard model and thus allow the decay branching ratio to lie in the 
range 2xl0~ 4 < B(b — > s+7) < 4.5 x 10 -4 . Another important aspect of mSUGRA is that 
it leads to the lightest neutralino to be the lightest supersymmetric particle (LSP) using 
the renormalization group analysis (see, e.g., Ref. ||33||) and hence with R parity invariance 
a candidate for cold dark matter. The current astrophysical data allows the amount of cold 
dark matter to lie in the range 0.1 < QcDAih 2 < 0.3[Q where flcDM = Pcdm/Pc where 
Pcdm is the matter density due to cold dark matter and p c is the critical matter density 
needed to close the universe, and h is the Hubble parameter in units of lOOkm/sMpc. In 
our analysis we assume the entire cold dark matter as arising from the lightest neutralino 
and thus impose the constraint on the neutralino relic density so that 0.1 < Q x oh 2 < 0.3. 

The result of the analysis is exhibited in Fig.l in the m — Mi plane, where m is the 
universal scalar mass and Mi is the universal gaugino mass, for fi positive and for values 
of tan/3 of 5,10,30 and 50. Thus yellow region-I at the top right hand side is disallowed 
either due to the absence of EWSB or due to the lighter chargino mass lying below the 
current experimental lower limit. Further, yellow region-II at the bottom is disallowed 
due to stau becoming the LSP. The red area on the boundary between the yellow and 
the white region is the zone satisfying the relic density constraint 0.1 < Q x oh 2 < 0.3 
(the blue filled circles are the additional allowed points if the relic density constraint 
includes the lower region 0.02 < Q x oh 2 < 0.1). As discussed above b — > s + 7 is not 
a strong constraint for positive [i and small tan/3. As a consequence the b — > s + 7 
constraint does not eliminate any parameter space for the cases tan /3 — 5 and tan j3 = 10 
as can be seen from Figs. (la) and (lb). However, for large values of tan/3 the b — > 
s + 7 constraint becomes important even for positive [i as can be seen from Figs.(lc) 
and (Id). We turn now to the imposition of the a e * p — a^ M constraint of Eq.(2). In 
the theoretical computation of the supersymmetric electroweak correction we take into 
account the full one loop supersymmetric contribution including the chargino-sneutrino 
and the neutralino-smuon exchanges. The result of imposing 2a variation around the 
central value gives the corridor shown in Figs.l by the upper (—2a) and the lower (+2<r) 
solid lines. Specifically the ranges of sparticle masses, Higgs boson masses and the [i 
parameter corresponding to la and 2a constraints are exhibited in Fig. 2. The spectra 
span a wide range and with the lower end of the scale corresponding to a low values of fine 
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tuning and the high end corresponding to high values of fine tuning (see, e.g., Ref. [ 
The spectrum of Fig. 2 leads us to the conclusion that sparticles satisfying la constraint 



are well within reach of the LHC|36| and some of the sparticles could also show up at 
RUNII of the Tevatron. The ranges allowed by the 2a constraint constraint have upper 
limits which in some cases reach as high as 3 TeV and may exceed the reach of the LHC. 
However, Fig. 2 does not contain the constraints of b — > s + 7 or of relic density. With 
these additional constraints the upper limits would be substantially smaller. 

We discuss now the direct detection of dark matter under the BNL2002 g M — 2 con- 
straint. The analysis of Fig.l shows that there are substantial regions of the parameter 
space where all the constraints, i.e., the g-2 constraint the 6^5 + 7 constraint and the 
relic density constraints are simultaneously satisfied. We discuss now the direct detection 
of dark matter under these constraints. There has been considerable work on the the- 
oretical predictions of the neutralino-proton cross section a xp which enters in the direct 



detection of dark matter and some recent literature can be found in Ref.[|37]]. Here we 
give the analysis of these cross sections in mSUGRA for positive /1 under the new g^ — 1 
constraint. Results are exhibited in Figure 3 for values of tan/? of 5,10,30 and 50. We 
find that the neutralino-proton cross sections consistent with all the constraints lie in the 
range 10~ 46 cm 2 < a xp < 2 x 10~ 43 cm 2 . A significant part of this range will be accessible 
at future dark matter detectors GENIUS KM and ZEPLINB9]. 



3 Conclusion 

The new Brookhaven measurement of — 2 has confirmed the earlier measurement with 
twice the precision. However, interpretation of what this implies for physics beyond the 
standard model is very sensitive to the errors in the standard model prediction. The main 
surprise at the end of last year was the switch in sign of the LbL contribution. There is now 
a general agreement that the sign of this contribution is positive and also an agreement 
between two independent evaluations |6|, 0] on its size. An important exception to this is 
the result of chiral perturbation theory, which although a more fundamental approach is 
beset by the appearance of unknown low energy constants which require a nonperturbative 
approach such as lattice gauge calculation for its evaluation. In the present analysis we 
have assumed the validity of the analysis on LbL of Knecht and Nyffeler and the revised 
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analysis of Hayakawa and Kinoshita. Taking account of all the errors the difference 
between experiment and the standard model is now predicted to lie between 1.6<r and 
2.6a. We have revisited the implications of the g M — 2 constraint for supersymmetry in 
view of the new result. Using the 2.Qa difference and a 2a error corridor our analysis points 
to significant regions in the mSUGRA parameter space consistent with the b — > s + 7, 
relic density and — 2 constraints. Within la error corridor all of the sparticles are 
accessible at the Large Hadron collider. Further, for low values of tan/? sparticles may 
also be accessible at Fermilab Tevatron. We also carried out an analysis of the neutralino- 
proton cross section er x <>_ p and find that a significant part of the parameter space would 
be explored by the future dark matter detectors such as CDMS (Soudan), GENIUS and 
ZEPLIN. Finally, the g^ — 2 experiment could lead to an even more stringent constraint 
if there was a reduction in the error associated with the standard model prediction. This 
could come about by eliminating the ambiguity on the LbL correction through a lattice 
gauge calculation analysis, and through improved low energy data on e + e~ — > hadrons 
in the analysis of hadronic vacuum polarization correction. 

We thank William Marciano for an interesting discussion. This work is supported in part 
by NSF grant PHY-9901057. 
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Figure Captions 

Figure 1: Analysis in mSUGRA exhibiting the allowed and the disallowed regions in the 
mo — Mi plane due to the radiative breaking of the electro- weak symmetry, the b — * s + 7 
constraint and the neutralino relic density constaint as discussed in the text for values 
of tan/3 of 5,10,30 and 50. The area to the right of the dotted blue line in each case 
is the allowed region by the b — > s + 7 constraint and the (red) dark area satisfy the 
relic density constraint 0.1 < Q x oh 2 < 0.3 (while the (blue) dark filled circles satisfy 
0.02 < Q x oh 2 < 0.1). The region I is discarded because of the absence of the electroweak 
radiative symmetry breaking or via the lower limit of chargino mass constraint. The re- 
gion II is eliminated because of stau becoming the LSP. The solid black lines exhibit the 
lower limit (a M (2o")) and the upper (a At (— 2cr)) limit corresponding to the ±2cr variation 
around the central value given by the first line of Eq.(2). The LEP 'higgs signal' is also 
indicated in the figures. 

Figure 2: Mass ranges of sparticles and Higgs bosons as well as the \i parameter in 
mSUGRA for fi > when the la and 2a — 2 constraints are imposed and tan (3 ranges 
from 2-60 and Aq = 0. The b — > s + 7 constraint and the relic density constraint are not 
imposed. For each (s)particle the vertical line to the left is the range in mass allowed by 
a 2a — 2 constraint while the vertical line to the right is the mass range allowed by a 
\a g^ — 2 constraint. 



Figure 3: A plot of the scalar proton-LSP cross-section a xp vs LSP mass m x which enters 
in the direct detection rates. The blue dots are the points allowed in mSUGRA and the 
red filled circles satisfy the constraint 0.1 < Vth 2 < 0.3. Black squares which form the 
black region satisfy the 2a constraints on b — > s + 7 and g^ — 2. The enclosed region on the 



top left hand side is the region where the DAMA collaboration^] claims a signal. The 



dot-dashed line on top is the upper limit from CDMS[41| while the solid line is the latest 
upper limit from EDELWEISS experiment [ VI \ . The lower dashed line is the sensitivity 
that will be reached by the CDMS experiment at Soudan mine while the dotted line is 



the sensitivity that will be achieved by the GENIUS detector[38 
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